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Abstract 

A general algorithm for handling the energy dependence of hadron-nucleon 
amplitudes in the nuclear medium, consistently with their density depen¬ 
dence, has been recently applied to antikaons, eta mesons and pions in¬ 
teracting with nuclei. Here we apply this approach to antiprotons below 
threshold, analyzing experimental results for antiprotonic atoms across the 
periodic table. It is also applied to antiproton and antineutron interactions 
with nuclei up to 400 MeV/c, comparing with elastic scattering and annihi¬ 
lation cross sections. The underlying pN scattering amplitudes are derived 
from the Paris NN potential, including in-medium modihcations. Emphasis 
is placed on the role of the P-wave amplitudes with respect to the repulsive 
S'-wave amplitudes. 

Keywords: antiproton-nucleon in-medium interaction, energy dependence, 
antiprotonic atoms, antiproton scattering, antiproton and antineutron 
annihilation, Paris NN potential 


1. Introduction 

The connection between hadron-nucleus empirical potentials near thresh¬ 
old and the underlying hadron-nucleon interactions has been studied for 
years by analyses of strong-interaction effects in exotic atoms and in stud¬ 
ies of elastic scattering of hadrons by nuclei n It was recognised in the 
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early 1970s that the 77-nucleus interaction near threshold was determined 
by the ^-nucleon scattering amplitude at subthreshold energies 0, H 0. 
Recently an algorithm was devised to account for the subthreshold energy 
dependence of the meson-nucleon amplitude in evaluating the meson-nucleus 
strong-interaction po tential in kaonic atoms and for strongly-bound K and 

Pionic atoms and elastic scattering 


r] mesons in nuclei 

of 22 MeV TT^ by nuclei have also been studied using this approach 

In the present work we apply this approach to the interaction of antipro¬ 
tons with nuclei near threshold. As in previous works, we are not interested 
in any single nuclear species but rather in global behavior. Therefore we 
handle only large data sets within ‘global’ comparisons between calculation 
and experiment, as was done for kaonic atoms [l^ and for pionic atoms 
and pion scattering 13|. Nevertheless, in order to assess the validity of the 


model, some annihilation cross sections are also considered. Due to the much 
stronger absorption of antiprotons in nuclei compared to pions, and even to 
antikaons, it is inevitable that ambiguities may exist in some of the conclu¬ 
sions. 

An extensive data-base for strong interaction effects in antiprotonic atoms 
is available from the PS209 collaboration at CERN 141. Results for elastic 


scattering of 48 MeV antiprotons on C, Ca and Pb nuclei are available from 


the pioneering experiments of the 1980s [l^ . For the free-space pN interac 


tion near threshold we used the 2009 version of the NN Paris potential 
This potential consists of a long-range one-pion and correlated two-pion ex¬ 
change terms, plus a short-range phenomenological term that includes an 
absorptive component representing pN annihilation. The potential param¬ 
eters are htted to some 4300 scattering data plus scattering lengths and 
scattering volumes extracted from antiprotonic hydrogen levels. This quali- 
hes the 2009 version of the NN Paris potential as a realistic potential. Other 
realistic, and ‘microscopic’ as well, NN potential models that have become 
available recently could, in principle, be used in p-nucleus calculations near 
threshold. These include (i) the Zhou-Timmermans model 0 which is also 
based on a long-range one-pion and correlated two-pion exchange terms, but 
uses a boundary condition description for its short-range term; and (ii) a 
Bonn-Jiilich NNLO chiral EFT potential model 18| with a similar long- 


range behavior that is subject, however, to a strict power counting hierarchy, 
and in which the short-range behavior is given by suitably determined con¬ 
tact terms. The present work is not intended to compare between different 
microscopic NN potential models, nor to study possible NN quasi-bound 
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states near threshold 


or 


^0 


[e.g. with quantum numbers ^^S'o |19 


or 


13 


Pn 20 


2l|h but rather to apply a microscopic model in the context of 


antiproton-nucleus interactions below threshold and at very low energies. As 
our handling of in-medium scattering amplitudes assumes some rather gen¬ 
eral properties of nuclei, we do not consider very light nuclei in the present 
work. 

NN potentials are found to be strongly attractive and absorptive in all 
the microscopic models known to us, including the various versions of the 
Paris potential and the recently published potentials mentioned above. This 
results, generally, in repulsive NN S'-wave scattering amplitudes at low ener¬ 
gies. Hence, the simple impulse-approximation optical potential, with 
tpAT the corresponding free-space pN t matrix and p the nuclear density, is 
repulsive at and near threshold. However, past global analyses of antipro¬ 
tonic atoms j^, ^ achieved good agreement with experiment when using 
an empirical (as opposed to ‘microscopic’) local attractive and absorptive op¬ 
tical potential, related to nuclear densities through a folded-in hnite-range 
interaction of rms radius about 1.1-1.2 fm. A strong density dependence 
of the effective, in-medium pN t matrix is apparently required in order to 
reverse the sign of the free-space pN t matrix in the medium and inflate its 
size, or a signihcant contribution from P-wave amplitudes is able to achieve 
it. This problem was recognised already in the 1980s, with several many- 


body mechanisms sug 
optical potential 12^ 


j;ested for obtaining an attractive low-energy p-nucleus 
but none of these works was able to test such 
proposed mechanisms in actual global analyses of antiprotonic atoms as the 
high-quality data of the PS209 experiment [l^ were non-existent then. Pre¬ 
vious attempts to add empirically a P-wave potential term or non linear 
density terms were unsuccessful, failing among other things to respect 
constraints imposed by neutron density distributions j^. More specihcally, 
an imaginary part of a P-wave potential compatible with an earlier version 
of the Paris potential could be accommodated, but then the real part 
of the S'-wave term was found to be incompatible with the Paris potential. 
It is therefore interesting to apply the present approach of treating energy 
and density dependence (reviewed in j29[) to the latest version of the Paris 
potential |16| . 

The paper is organized as follows. Section [2] is a brief description of the 
present approach. In Subsection 12. II we introduce the in-medium kinematics 
satished in hadron-nucleus collisions, in which the hadron-nucleon center of 
mass (cm) energy depends also on the momenta of the participating par- 
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tides, and in Subsection 12.21 we discuss in-medium corrections to the free 
pN amplitudes such as Pauli correlations. Section [3] describes the free-space 
S'-wave and P-wave pN input amplitudes derived from the 2009 version of 
the NN Paris potential and used in the present work. Low-energy Np an¬ 
nihilation is also considered in this section. Section 0] reports on results of 
comprehensive hts to strong-interaction observables in antiprotonic atoms, 
with special emphasis placed on the role played by the P-wave pN ampli¬ 
tudes in reproducing the main features of best-fit empirical p-nucleus optical 
potentials. In Section |5] we present some results for the elastic scattering of 
48 MeV antiprotons by C, Ca and Pb nuclei. Annihilation cross sections of 
antiprotons and antineutrons on nuclei are also briefly mentioned. Section [3] 
offers a brief discussion and summary of the present study. 


2. Theoretical background 

Here we present the essentials of the theoretical background for the present 
work, referring to previous publications for further details. 

Strong interaction observables in pionic and kaonic atoms are usually 
calculated from the relativistic Klein-Gordon equation 

- 2/i(P + Kpt + 14) + (14 + Bf] = 0, (1) 


where h = c = 1 is implicitly assumed. Here p is the meson-nucleus reduced 
mass, B is the complex binding energy, 14 is the hnite-size Coulomb inter¬ 
action of the meson with the nucleus, including vacuum-polarization terms. 
14 pt is the optical potential describing the strong interaction of the meson 
with the nucleus. For antiprotonic atoms we can use this equation because for 
a given I in good approximation it gives the spin-averaged results of the Dirac 


equation j22|. This is certainly acceptable with the experimental results of 


the PS209 collaboration which forms the basis of the present work. 

One of the aims of the present work is to elucidate the role of the P-wave 
part of the anti pro ton-nucleon interaction in antiproton-nucleus interactions 
near threshold [271, |28|, Therefore we include it explicitly in the optical 
potential, which we take in analogy to the pion-nucleus potential 


2ppVopt{r) = q{r) + 3V • a(r)V, 


( 2 ) 


where pp is the p-nucleus reduced mass but unlike for pions, a factor 2/-I-1 = 3 
is introduced explicitly into the P-wave part to match the normalization of 
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the amplitudes of the following Section. The S'-wave part is written as 

g(r) = -47r(l + —){bo[pn{r) + Pp{r)] + bi[pn{r) - Pp{r)]} (3) 

rriN 

and the P-wave part 

a{r) = 47r(l + —)~\co[pn{r) + pp{r)] + Ci[pn{r) - Pp(r)]}, (4) 

rriN 

where p„(r) and Pp{r) are the neutron and proton densities normalized to 
N and Z, respectively, with N+Z=A. If evaluated at threshold the parame¬ 
ters 6o,i and co,i are related to the scattering lengths and to the scattering 
volnmes, respectively. However, in the present work these parameters are 
evalnated at density-dependent energies as explained below. 

2.1. In-medium kinematics 

The model underlying the snbthreshold energy algorithm adopts the Man¬ 
delstam variable s = {Ep + — {pp + as the argument transform¬ 

ing free-space to in-medium antiproton-nucleon amplitndes, where both the 
p and the nucleon variables are determined independently by the respective 
environment of a p atom and a nncleus. Conseqnently, unlike in the two-body 
cm system, here pp pn does not vanish, and one gets to a good approxi¬ 
mation (pp -|- p]\fY Pp + P% upon averaging over angles. The energies are 
given by 

Ep = mp- B^, Epf = mN- Bn, (5) 

where B^ is the real part of the p binding energy in the atom. Bn the (real) 
binding energy of the nncleon and m are masses. For the p momentnm we 
snbstitnte locally 

2 

^ = -B^ - Re Ropt - Re (6) 

^11 Lp 

For the nucleon we adopt the Fermi gas model (FGM), yielding in the local 
density approximation 

( 7 ) 

where p is the local density, p is the average nnclear density and Tn is the 
average nncleon kinetic energy which assnmes the valne 23 MeV in the FGM. 
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Defining S^/s = ^/s — with = mp + m^v, then to hrst order in 
B/Eth and (p/Eth)^ one gets 

= -Br,p/p - UTn{p/pY^^ + B^p/po] + Kpt + K(p/po)'/'], (8) 

with ^jv = ^N/i'iTT'N + ^p), ^p = mp/{mM + mp), and po = 0-17 fm”^. 
Following previous applications 0,0 H B113 an average binding energy 
value of -Btv = 8.5 MeV is used. The specific p/Pn and p/p forms of density 
dependence ensure that 5\fs —)• 0 when p —)■ 0 lOl . 

Another variant of Eq. (jS]) is obtained when considering the minimal 
substitution requirement, the importance of which for incorporating electro¬ 
magnetism in a gauge-invariant way into the pion optical potential was hrst 
pointed out by Ericson and Tauscher Ml] and more recently emphasized by 
Kolomeitsev, Kaiser and Weise . Indeed, the application of minimal sub¬ 
stitution has been successful in analyses of pionic atoms 0, 33] and pion 


Here E = Ep + Ejq is replaced by E — 14 


scattering at low energies [13|, [34 
and then Eq. ([8]) becomes 

= -Bnp/p - iN[TN{p/pf'^ + B^p/po + Kc(p/po)'/'] + ^pRe Ropt- (9) 

Eq. (]9]) is used in the present work to handle the in-medium kinematics in 
antiprotonic atoms. At energies above threshold, in applications to scattering 
and in-hight annihilation, the term —^pfBpp/ po is replaced by ^ArEiab, with 
Eiab the beam kinetic energy, leading to 

hVi = -Bnp/p - ^n[Tn{p/p)^^^ - Eiab + Rc(p/po)'/'] + epRe Ropt- (10) 
2.2. In-medium amplitudes 

To obtain in-medium amplitudes from the free space ones we apply the 
multiple scattering approach of Waas, Rho and Weise (WRW) as used 
for kaonic atoms in Ref. [13 . Since the isospin structure of KN and pN is 
the same, we use for the R-wave potential q{r) of Eq. ([3]) the form given in 
Ref. [13, 


g(r) = —dvr 


(2/pp fpn) 2pP 


+ 


fpn ( 2 Pp R 


1 R ji^k=oi/^fpp fpn)pij’') 1 R ^^k=ofpnPij'') 


( 11 ) 


where p = pp + pn and f = (f is a. p-nucleus cm amplitude related to the 
two-body cm amplitude /, with ( given by 

A 1 Pp 


C = 1R 


A 


mjv 


( 12 ) 
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To leading order, ^^=0 = 'with pp the local Fermi momentum, account¬ 

ing for nuclear Pauli correlations. These WRW medium corrections apply 
only to the S'-wave part of the potential j^. The P-wave part, Eq. (|1]), 
is taken from the free amplitudes at energies prescribed by Eq. (|9]) for an¬ 
tiprotonic atom applications, or by Eq. ffTOj) for applications above threshold, 
without incorporating further nuclear correlations such as realized in pionic 
atoms by the Lorentz-Lorenz modihcations. This is justihed in p atoms owing 
to the extremely low nuclear densities encountered in their analysis. 



P/Po 


Figure 1: Examples of the density to energy transformation from Eqs. and (nni), using 
Eq. ([TT]). 


Equations ([9]) and (fTOj) dehne a density to energy transformation through 
the scattering amplitudes used in the calculation of antiproton-nucleus op¬ 
tical potential Eq. (nj. As the real part of the potential determines the 
energy and, in turn, the energy and density determine the amplitudes, a self- 
consistent solution is required. Examples for this transformation are shown 
in Fig. [1] for p atoms of ^°Ca and for elastic scattering of 48 MeV p using 
the amplitudes from the 2009 version of the Paris NN potential. The energy 
shift 5y/s for p atoms comes out negative dehnite over the full density range 
considered, reaching fairly large values of up to about —60 MeV at nuclear- 
matter density po) iii agreement with its behavior for kaonic atoms [l^. For 


7 




low-energy scattering, in contrast, 6y/s changes sign as a function of the 
density from positive values to negative ones around 0.6po, and its slope ex¬ 
hibits a marked discontinuity caused by the in-medium S'-wave amplitude, 
Eq. flTT|) . switching from repulsion to attraction as the density is increased. 
In the case of p atoms, where the input pN S'-wave amplitudes are apprecia¬ 
bly larger in size than for 48 MeV antiprotons, the transition from repulsion 
to attraction takes place at very small densities, thus making it unobservable 
on the curve. Further discussion of the WRW mechanism responsible 
for the in-medium sign reversal of the real part of the S-wave potential ^(r), 
Eq. (ITT]) , is deferred to Appendix A, including some realistic estimates. 

3. Free-space input amplitudes 

3.1. General 



Figure 2: Antiproton-nucleon S'-wave scattering amplitudes around threshold from the 
2009 Paris potential 0. Negative values represent attraction and absorption. 


Here we briefly describe the free-space pN amplitudes derived from the 
2009 version of the NN Paris potential [l^. It was suggested by Green and 
Wycech that the pN amplitudes most appropriate for use in antipro¬ 
tonic atoms are half off-shell, (p \f{E)\p') with one on-shell cm momentum 







p = y/mis[E and one off-shell cm momentnm which we choose as p' = 0 to 
focus on near-threshold energies. Figs. |2] and |3] show such F-wave and P-wave 
half off-shell pN amplitudes, respectively, derived from the 2009 version of 
the NN Paris potential. The pn amplitudes are pure isospin T = 1 NN 
amplitudes, whereas the pp amplitudes are equal-weight mixtures of isospins 
T = 0 and T = 1 iVA^ amplitudes. Apart from their isospin structure, these 
amplitudes represent angular-momentum averages, appropriate for use in an¬ 
tiprotonic atoms, over states denoted by ( 2 'r+i)( 2 S'+i)p^^ where J, S, L are the 
total, spin and orbital angular momentum, respectively. The actual pp and 
pn A-wave (L = 0) and P-wave {L = 1) amplitudes plotted in Figs. |2]and|3l 
respectively, are derived from the T = 0 and T = 1 pure isospin ampli¬ 
tudes, which are obtained by the appropriate angular-momentum averaging 
of £xed-T amplitudes as follows 


ft 


Esj +1) fli 

Esj (2^ +1) ' 


(13) 



Figure 3: Antiproton-nucleon P-wave scattering amplitudes around threshold from the 
2009 Paris potential [l^ . Negative values represent attraction and absorption 


The real part of the Paris potential is predominantly attractive and gen¬ 
erates quasi-bound states or resonances in specihc partial waves. However, 
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in the S'-wave amplitudes these are washed out by the very strong absorp¬ 
tion which acts repulsively. In particular the 2009 version of this potential 
generates a broad quasi-bound ^^S'o state a.t E = —4.8 — *26 MeV. This 
state has a small statistical weight and it causes a small anomaly just below 
threshold, producing no effect on ^nuclear observables near threshold. It 
is seen from Fig. [2] that the free S'-wave amplitudes, in distinction from the 
underlying attractive potential, represent repulsion over a very wide energy 
range. This disagrees with the empirical attractive potential deduced from 
previous analyses of antiprotonic atoms unless medium effects, such as the 
WRW prescription discussed in Section |2] or other amplitudes, e.g. P-wave, 
produce attraction. 

Fig. in] shows the P-wave amplitudes from the Paris potential. In the 
P-wave states the short range annihilation has much weaker effect and the 
attractive nature of the potential is evident. In particular a quasi-bound, 
fairly narrow state at P = —4.5 — *9.0 MeV is generated in the ^^Pi partial 
wave. Due to its sizable statistical weight it is seen clearly in the hgure. 
In contrast to the S'-wave quasi-bound state discussed above the P level is 
robust, arising in all versions of the Paris potential. In the 2009 version its 
position is rather reliably hxed by the pp scattering volumes extracted from 
the antiprotonic hydrogen atom. 


3.2. Low energy antinucleon-proton annihilation 

Before applying the Paris 2009 amplitudes to the p-nucleus interaction 
close to threshold we examine how well these amplitudes describe experi¬ 
mental results for the iV-proton interaction at very low energies. 

The obvious starting point is the strong interaction shift and width ob¬ 
served in pH atoms 36l. l37| which, with the help of the Deser 
equations 


38 


39 


or Trueman 

approximates the scattering length for the pp interaction. Us¬ 
ing a modern version of the latter ^ we have calculated also the scattering 
length in the absence of the Coulomb interaction using Eq. (20) of Ref. j4lj |. 
The result of —0.81-l-*0.72 fm is in reasonable agreement with the amplitudes 
used in the present work. 

Next we compare predictions of in-flight annihilation cross sections up 
to 400 MeV/c with calculations based on the S- and P-wave amplitudes. 
Figure m shows np annihilation cross sections from the OBELIX collaboration 


42j . As some of the quoted experimental errors are unrealistically small, we 


have assumed errors of ±10% at all points. No parameter adjustment has 
been made in the calculation. Except for the small bump near 200 MeV/c the 
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(MeV/c) 


Figure 4: Comparison between predictions and experiment 4^ for np annihilation 



Figure 5: Comparison between predictions and experiment 
lation. 


3, S ES El for pp annihi- 
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agreement is very satisfactory. Fig. Elshows similar results for pp annihilation 


431144 l45l. |46| . A subtle point may arise here because the 2009 amplitudes 


were calculated without the Coulomb interaction. An approximate correction 
has therefore been applied by multiplying the F-wave and P-wave amplitudes 
by the Coulomb phase correction, namely, and , respectively, where 
a are the Coulomb phases. At the lowest momentum the calculated cross 
section is then increased by 28%, whereas near 400 MeV/c this correction is 
less than 2%. The overall agreement between predictions and experiment is 
certainly acceptable. 


4. Antiprotonic atoms 

The PS209 data used in the present work consist of strong interaction 
observables for 27 nuclear targets from “^^Ca to totaling 84 data points 

(l^ . Six points from earlier measurements on were also included. 

Typical examples of results for various choices of potentials are shown in 
Table [H Values of for fho 90 points indicate the quality of £t, obtained 
when varying 2-3 parameters. Wher e appro priate, the present results are in 


full agreement with earlier analyses j22|, |23 


The hrst row of Table [T] shows, as a reference, the best £t obtained with 
a local emp irical attractive and absorptive potential, including hnite-range 


folding [22|, l23[ . The very strong absorption of antiprotons in nuclei conhnes 


the interaction with the atomic p to the surface of the nucleus and this is 
demonstrated in the second row of the table, where we used only an empirical, 
energy independent P-wave potential. It is evident that equally good fits can 
be obtained when the potential is centered near the surface, where the gradi¬ 
ent terms are effective. The third row is for a potential constructed from the 
free-space Paris 2009 amplitudes of the previous section, in which the WRW 
in-medium modification flTT]) is incorporated with in-medium kinematics sat¬ 
isfying the 5y/s algorithm for p atoms, Eq. ([9]). No adjustable parameters are 
included. It is clear from the value of fh^t no agreement with experiment 
is possible. The reasons for this failure are further discussed below in the 
last paragraph of the present section. 

The fourth row of Table [Ushows the results obtained by replacing the free- 
space P-wave term by an empirical one, while retaining the F-wave part as 
above, namely, the medium-modihed amplitudes at density-dependent ener¬ 
gies given in the present model. The fit to the data, varying two parameters, 
is as good as the best one. Comparing parameters with the results on the 
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Table 1: Comparisons between calculation and experiment for p atoms, using various 
options of S'-wave and P-wave potentials. The symbol ’09 means that pN amplitudes 
were used over a range of ^/s values as given by Eq. ([9|) , including the WRW modification 
Eq. (Ilip . When indicated, these were multiplied by a scaling factor given in parentheses. 
When units of fm^ are listed, the corresponding parameter was empirical. See text for 
more details. 



S'-wave 

Real P-wave 

Imag. P-wave 

X^(90) 

1 

emp. 

- 

- 

199 

2 

- 

1.9T0.1 fm^ 

2.8±0.1 fm^ 

206 

3 

’09 

’09 

’09 

2304 

4 

’09 

2.9±0.1 fm^ 

1.8±0.1 fm^ 

203 

5 

’09 

’09x(-10.0±0.9) 

’09x(3.1±0.2) 

571 

6 

’09 

2.9T0.1 fm^ 

’09x(1.3±0.1) 

218 


second row, we see that the real P-wave part is more attractive now as it has 
to overcome the repulsion of the S'-wave part near the surface. The opposite 
is true for the imaginary part because the imaginary part of the S'-wave term 
already provides some of the absorption. 

Row hve of the table shows an attempt to modify the free-space P-wave 
amplitudes by a factor, separately for the real and imaginary potentials, 
within the present full approach. This was not successful as is seen from the 
value of ^ind by the scaling factor imposed on the real part which requires 
a sign change and an order of magnitude enhancement. Realizing that the 
major difficulty with the free-space amplitudes is probably with the real part 
of the P-wave, the last row shows a very good £t to the data when rescaling 
the imaginary part of the P-wave amplitude by 30% while using an empirical 
real part. The latter agrees well with the one in row 4. 

Figures [6] and [7] show real and imaginary parts of the p-^°Ca potentials, 
respectively, typical of the results for the whole data base. Solid curves show 
the best-£t empirical potential which is attractive and obviously absorptive. 
The repulsive real part obtained from the free-space amplitudes at thresh¬ 
old is also shown (dot-dash) and the effects of the medium modihcation are 
clearly seen on the dashed curves. In the medium the WRW modihcation, 
assisted by the negative 5y/s shift for p atoms (see Fig. [1]), turns the repulsion 
into attraction as demonstrated in Appendix A, but far down the surface the 
effect disappears and the potential becomes again repulsive, differing signif- 
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r(fm) 

Figure 6: Real part of antiproton-nucleus potentials for ^°Ca p atoms. Solid curve for 
the best-fit empirical potential, dot-dash curve for the unmodified S-wave potential at 
threshold, dashed curve for the medium-modified S-wave potential. Vertical dotted line 
marks the half-density radius of ^°Ca. 



Figure 7: Imaginary part of antiproton-nucleus potentials for ^°Ca p atoms, see caption 
to Fig. [6] for the meaning of the various curves. 
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icantly from the empirical potential that hts well the data. Regarding the 
in-medium imaginary potential, it falls short of its empirical counterpart in 
the surface region. The WRW medium modihcation reduces the imaginary 
potential in the interior by a factor about three while enhancing it in the 
surface, but insufficiently to match the best-£t potential. The vertical dot¬ 
ted line in both hgures indicates the half-density radius of the nucleus, thus 
demonstrating that the relevant region is at larger radii. With these observa¬ 
tions it is clear why the application of the medium-modihed Paris potential 
amplitudes to antiprotonic atoms fails, as was shown in the third row of 
Table [T] 


5. Af-nucleus above threshold 


Calculated in-flight annihilation cross sections for p on light muclei, using 
the present in-medium procedures, are in good agreement with the very few 


available experimental results 1^ . However, the best test of a potential 


is by comparing predictions to measured differential cross sections. Although 
at 48 MeV beam energy S and P waves might be insufficient to fully describe 
microscopically the interaction of p with medium weight and heavy nuclei, 
we report below some features observed using the Paris 2009 potential at 
that energy. 

5.1. Scattering of 48 MeV antiprotons 

Measurements of elastic scattering of antiprotons by ^°Ca and ^osp]^ 
were made in the 1980s and analyzed using standard low-energy optical model 
methods, see Janouin et ah 1^ and references therein. Here we apply the 


approach used above for antiprotonic atoms also to the p-nucleus interaction 
at 48 MeV energy. In parallel with p atoms, excellent £t to the data could be 
obtained with an empirical local potential, including hnite-range folding with 
rms radius of 1.5±0.1 fm, as is summarized in the hrst row of Table |2l This 
range is somewhat larger than the corresponding one for p atoms. An attempt 
to use only an empirical P-wave potential was not successful unless hnite- 
range folding was introduced, unlike with atoms, and then only a moderate £t 
was obtained, see the second row of the table. Very good hts were possible by 
a combination of empirical local and P-wave potentials, but then correlations 
prevented achieving a unique solution. Experience shows that no more than 
three meaningful parameters can be derived here. 
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Table 2: Comparisons between calculation and experiment for 48 MeV p scattering by 
12c, 4°Ca and ^ospb, using various options of S'-wave and P-wave potentials. The symbol 
’09 means that pN amplitudes were used over a range of ^/s values as given by Eq. (fTOll . 
including the WRW modification Eq. (HB). When indicated, these were multiplied by 
a scaling factor given in parentheses. When units of fm^ are listed, the corresponding 
parameter was empirical. The last column lists the rms radii of the finite-range folding 
applied to the adjusted terms. See text for more details. 



P-wave 

Real P-wave 

Imag. P-wave 

X^(83) 

range (fm) 

1 

emp. 

- 

- 

151 

1.5±0.1 

2 

- 

0.5±0.05 fm^ 

0.8±0.05 fm^ 

201 

0.9±0.1 

3 

’09 

’09 

’09 

1823 

- 

4 

’09 

0.31±0.03 fm^ 

0.50T0.03 fm^ 

176 

1.5±0.1 

5 

’09 

’09x(-2.1±0.2) 

’09x(2.4±0.2) 

495 

- 

6 

’09 

’09x(-1.92±0.04) 

’09x(1.32±0.09) 

218 

1.5±0.1 


The third row of Table [2] is for using in-medium amplitudes generated 
from free-space amplitudes of the Paris potential by using in-medium kine¬ 
matics, Eq. m, and applying the WRW modihcation flTT]) to the S-wave 
part without any adjustable parameters. Clearly that is unacceptable, but 
row 4 shows a good £t when the 2009 energy dependent Paris-potential P- 
wave amplitudes are replaced by an empirical P-wave term. Row 5 is for 
attempts to £t the data by applying scaling factors to the Paris potential 
P-wave amplitudes. No £t is possible but we note that the real part required 
a change of sign, as is the case with atoms (see Table [T]). Finally, row 6 
shows that an almost acceptable £t is possible with hnite range applied to 
the P-wave term but with signihcant scaling factors. 

Figure [8] shows, as an example, the imaginary part of p nuclear potentials 
for 48 MeV antiprotons interacting with ^‘^Ca. As for p atoms, there are sig¬ 
nihcant differences in the surface region between the in-medium microscopic 
P-wave potential and the best-fit empirical one. Figure [9] compares experi¬ 
mental differential cross sections for elastic scattering of 48 MeV p from ^°Ca 
with calculations made with the best-fit empirical potential (row 1 of Ta¬ 
ble [2]) and with the in-medium microscopic potential obtained from the Paris 
2009 free-space amplitudes (row 3 of the table). The poor agreement near 
the hrst minimum indicates that the microscopic potential is inadequate near 
the surface region of the nucleus. Indeed it is seen from Fig. [3] that whereas it 
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Figure 8: Imaginary part of antiproton-nucleus potentials for ^°Ca at 48 MeV. Solid curve 
for the best-fit empirical potential, dot-dash curve for the unmodified S'-wave potential 
at 48 MeV, dashed curve for the medium-modified S'-wave potential. Vertical dotted line 
marks the half-density radius of "^^Ca. 



Sam. (deg.) 


Figure 9: Differential cross sections for elastic scattering of 47.8 MeV p by ^°Ca 15j. Solid 
curve for the best-fit empirical potential (row 1 of Table dashed curve for 2009 Paris 
potential (row 3 of Table [2]). 
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might be possible for the P-wave part to add the extra absorption required, 
there is no chance that the real part of the P-wave amplitude will close the 
gap between empirical potentials and the in-medium microscopic potentials 
generated from the free-space amplitudes of the Paris potential. That is 
also the conclusion reached from inspection of Table |2l However, it must be 
emphasized that we have retained all along the S'-wave amplitudes as given 
by the Paris potential. Therefore the difficulties with the P-wave part es¬ 
sentially mean that there are some inconsistencies between the two types of 
amplitudes, within the present model of handling amplitudes in the nuclear 
medium. A need to include a D wave above threshold in a microscopic model 
cannot be ruled out. 


5.2. Annihilation on nuclei 

There has been only a handful of measurements of annihilation cross sec¬ 
tions of antiprotons on medium-weight and heavy nuclei at energies close to 
threshold, where the present optical model approach is being tested. In con¬ 
trast, measurements of annihilation cross sections of antineutrons on nuclei 
across the periodic table at seven momenta between 76 and 375 MeV/c were 


made by Astrua et ah [50[. These results have been compared [23|, 1^ with 


predictions by empirical optical potentials that £t quite well the few avail¬ 
able p-nucleus annihilation cross sections. Here we compare these results 
with predictions by potentials based on the Paris 2009 amplitudes. 

Figure [TUI shows, as an example, experimental annihilation cross sections 
for antineutrons on Sn (solid circles) and a single point (open square) for 
annihilation of antiprotons on the same target |52|. Calculations are shown 
as solid and dashed lines for the full 6y/s model and for 6y/s = .^Ari^iab, 
respectively. The weak sensitivity to the model is due to the interaction being 
conhned to the extreme surface region of the nucleus. The agreement with 
experiment for the single p point is very good. The sharp disagreement for the 


antineutrons was dicussed in j23|, , based on empirical optical potentials. 


It persists also here when the potential is constructed from the Paris 2009 
amplitudes. Very recently it was suggested by Bianconi et al. j^, using 
qualitative arguments, that annihilation cross sections could be enhanced 
by the interaction between the Coulomb held of a nucleus and an induced 
electric dipole of the antineutron. We have tested this idea using a realistic 
polarizability of 0.001 fm^ for the antineutron but the effect is totally 
negligible. Increasing the polarizability by an order of magnitude will lead 
to effects smaller than 1%. 
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Figure 10: Comparisons between calculated and experimental annihilation cross sections 
on Sn. Solid circles for n , open square for p ■ Solid lines are for using in-medium 
kinematics, Eq. m, and applying the WRW modification Eq. (HU. Dashed lines are for 
= ^TvElab- 


6. Discussion and summary 

The ability of a simple empirical optical model approach to describe well 
the interaction of sub-threshold and of low energy antiprotons with nuclei has 
been known for some time. Among other things, this success could result from 
the interaction being confined to the extreme surface region of the nucleus 
due to the very strong absorption of p and n in nuclear matter. Nevertheless, 
it was interesting to see how well can more microscopic approaches do in this 
respect. This has been done in the present work using the 2009 Paris NN 
potential fl^ . 

Inspection of optical potentials based on a recent algorithm for handling 
microscopic scattering amplitudes in the nuclear medium shows that it is nec¬ 
essary to include contributions from the P-wave part of the pN interaction, 
as suggested previously [^. In the present case this is necessary because 
the S'-wave amplitudes are repulsive throughout the full energy range, as 
seen in Fig. |2l and although the WRW medium modification makes them 
attractive in the interior, they are nevertheless repulsive at the relevant low 
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density region of the nnclens. This applies both to p atoms and to scattering 
at 48 MeV beam energy. The reversal of sign of the real potential in a snffi- 
ciently dense matter is exclnsively dne to the in-medinm WRW modihcation 
flTT]) and not dne to the in-medinm kinematics algorithm, Eqs. (jO]) and ffTOj) . 
Tests show that this phenomenon is linked to the particularly large values 
of both real and imaginary part of the S'-wave amplitudes and it disappears 
when the amplitudes are reduced by a factor 3 or more. In practice we hnd 
that the real part of the P-wave amplitudes cannot substitute for the missing 
S'-wave attraction without further sizable modifications, presumably because 
they change sign sharply very close to threshold (Fig. [3]) which is the relevant 
energy in the present studies. 

We have also studied possible shifts, Eq. ([2]), in the energy where the 
P-wave amplitudes are evaluated and found, in fact, that shifting to about 
8 MeV lower energies greatly improves the agreement with experiment. How¬ 
ever, judging by the pp and np annihilation cross sections. Figs. H] and [5] 
respectively, the energy scale cannot be changed by more than ±1 MeV. 

Finally, the large discrepancies between measured and calculated fi-nucleus 


annihilation cross sections, observed earlier with empirical potentials [23|, iM 
persist in the present work based on a more microscopic approach. 
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Appendix A. Prom repulsive free-space pN amplitudes to attrac¬ 
tive in-medium pN amplitudes 

Here we demonstrate the WRW mechanism, Eq. OB. responsible for 
turning repulsive free-space pN S'-wave amplitudes into attractive, density 
dependent, in-medium pN amplitudes. It is instructive to estimate the size 
of the free-space pN repulsive S'-wave cm amplitudes / required to make the 
corresponding in-medium amplitudes /medium(p) attractive somewhere at the 
nuclear surface and inward, say beginning at density p = po/8. For simplicity, 
we assume that fpp = fpn, so that Eq. flTT]) at this particular value of density 
assumes the form 

/ 

/medium(p = Po/8) = \ 0 Q5 f ^ (^•^) 
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where / = /j? + ifi is given in fm. If / is purely real (// = 0) and re¬ 
pulsive (/ij < 0), then /medium(p > Po/8) > 0, Corresponding to in-medium 
attraction, for a sufficiently strong repulsive free-space amplitude with val¬ 
ues of fn < -1.54 fm. If / is purely absorptive {/r = 0), then /medium(p) 
is attractive and absorptive everywhere for any nonzero value of //. Tak¬ 
ing a representative subthreshold value of // = 1 fm from Fig. [21 one 
gets /medmm(p = Po/8) = (0.46 -|- *0.70) fm corresponding to a medium- 
size attraction and a somewhat reduced absorptivity. Increasing fj, the 
in-medium attraction at po/8 increases up to a maximum value and the 
in-medium absorptivity continues to decrease. Finally, if / is equally re¬ 
pulsive and absorptive {—/r = fi), Re/medmm(p > Po/8) > 0 for values of 
Jr < —0.77 fm. Taking again a representative value of // = 1 fm, one gets 
/medium(p = Po/8) = (0.55*1.83) fm, again corresponding to a medium-size 
attraction, but to a sizable absorptivity enhanced by almost a factor of 2. 
Note from Fig. [2]that in order to satisfy the minimum-repulsion requirement 
of Ir < —0.77 fm, one has to go slightly below threshold. 
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